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We probe the three dimensional flow structure and rheology of gravitational (non-density matched) 
suspensions for a range of driving rates in a split-bottom geometry. We establish that for sufficiently 
slow flows, the suspension flows as if it were a dry granular medium, and conflrm recent theoretical 
modelling on the rheology of spht-bottom flows. For faster driving, the flow behavior is shown to 
be consistent with the rheological behavior predicted by the recently developed "inertial number" 
approaches for suspension flows. 
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Flows of granular materials submersed in a liquid of 
unequal density have started to attract considerable at- 
tention [IH5] and are relevant in many practical appli- 
cations [5]. These materials, which we will refer to as 
"gravitational" suspensions, clearly differ from density 
matched suspensions, which have been studied in great 
detail [THIO] . Gravitational suspensions exhibit sedimen- 
tation, large packing fractions and jamming of the mate- 
rial, which suggests a description similar to dry granular 
matter HUES]. 

In the last two decades, various flow regimes have been 
identified for dry granular matter. Sufficiently slow flows 
are frictional: the ratio of shear (driving) to normal (con- 
fining) stresses becomes independent of flow rate if the 
material is allowed to dilate [12l [E]. Faster flows are 
referred to as inertial: here the effective friction coef- 
ficient depends on the so-called "inertial" number /, 
which is a non-dimensional measure of the local flow 
rate [Hill [15]. 

For gravitational suspensions, the presence of liquid 
instead of gas as interstitial medium strongly affects the 
microscopic picture — how should we think of the flow 
of such suspensions? Pouliquen and coworkers proposed 
that the ratio of the strain rate and settling time. Is, 
would play a similar role as the inertial number in dry 
granular flows [5,. They furthermore conjectured a de- 
pendence of the effective friction coefficient fi on Is sim- 
ilar to the dry case, and applied this rheological law to 
capture the behavior of underwater avalanches jl6| . 

Here we test this picture by combining 3D imaging and 
rheological measurements of the flow of gravitational sus- 
pensions in a so-called split-bottom geometry (Fig. 1). 
This geometry has two main advantages. First, the flow 
rate, which is the key control parameter in the inertial 
number framework, can be varied over several orders of 
magnitude, allowing us to access slow flows as seen in 
plane shear |3l [17] , faster flows as seen in gravity driven 
flows US] , and the crossover regime in between - some- 
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FIG. 1: (Color online) (a) The experimental setup used for 
flow visualization in a '"split-bottom"' geometry, (b) Exam- 
ple image of a single cross section, displaying half the box. 
(c) Geometry used for rheological measurements. (d) Geome- 
try used for flow visualization. 



thing not achieved in previous studies of gravitational 
suspensions (21 [T71 [TS] . Second, extensive experimental 
and numerical work [TUH21] has shown that the split bot- 
tom geometry produces highly nontrivial slow dry gran- 
ular flows. A simple frictional picture is not sufficient to 
capture these flows p51l26j . so that testing whether these 
proflles also arise in slowly sheared gravitational suspen- 
sions is a stringent test for similarities between slow dry 
flows and slow gravitational suspension flows. 

Setup — The split-bottom geometry is sketched in 
Fig. [T^, and consists of a square box, 15 cm in width 
with transparent acrylic walls, at the (rough) bottom of 
which a (rough) disk of radius Rg = 4.5 cm can be ro- 
tated at rate ft. 

We use monodisperse acrylic spheres with a diameter 
d of 4.6 mm (Engineering Labs); all our results are qual- 
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itatively similar for 3.2 mm particles. The particles are 
suspended in a mixture of some 78% Triton X-100, 13 % 
water, and 9% ZnC12 (by weight) [31] with a fluorescent 
dye added (Nile Blue 690) . The refractive indices of par- 
ticles and fluid are approximately 1.49 and match closely 
— we adapted the recipe from [2H]- The fluids viscos- 
ity rjf is 0.3 (±0.05) Pa s, and the difference in density 
between the fluid and the particles is about 100 kg/m"^. 

The particle motion is visualized by illuminating the 
suspension with a thin ( < 200/im ) laser sheet [311291 130]. 
The laser (Stocker Yale, 635 nm) is aligned parallel to the 
bottom of the box (Fig. [T^) and mounted on a z-stage 
which allows the illumination of slices of the suspension 
at different heights h [30] . Image acquisition is done with 
a triggered 12 bit cooled CCD camera, and contrast is 
sufficient to image half of the box (Fig. ^p) . We use a 
Particle Image Velocimetry-like (PIV) method to obtain 
the normalized azimuthal velocity w(r, z) = ve/{rfl) in 
slices of constant z. Combining these slices, we recon- 
struct the full angular velocity field as function of radius 
and depth for a range of driving rates. An overview of 
the imaging technique will be published elsewhere [31) . 

Rheological experiments were carried out by driving 
the disk from above with a rheometer (Anton Paar MCR 
501) — see Fig.jTj;. Velocimetry measurements were done 
by driving the disk from below with a DC motor — see 
Fig. [iji. There is always at least half a centimeter of fluid 
above the suspension to ensure that the surface tension 
of the fluid will not affect the dilation of the packing. 

Constitutive equation — We derive the constitutive 
equation for our suspension in the modified "iner- 
tial number" approach [2]. The typical rearrangement 
timescale for the particles in the suspension, given the 
viscosity and relative density of the particles becomes 
d/vini — Vf/Pctj where v-mf, P and a are settling ve- 
locity, pressure and porosity, so that the 'inertial num- 
ber' becomes: Is = -p^ [1]. The shear stress stress r 
is then written as r = iJ,{I)P^ with fi{I) an empirical 
friction function. For small /, fi{I) can be expanded: 
^{I) = ^0 + fill [2 [32], with /iQ and /ii empirical 
values. Combining this with the expression for Is, we 
arrive at: 




FIG. 2; (Color online) Frictional to viscous crossover in flow 
profiles, (a) The predicted flow field for dry granular flows 
!jjD{r,z) from Refs. |19H25j . (b-e) Measured velocity fields 
uis{r,z) at driving rates f2 = 8.3 x 10~^ rps (b), D, = 8.3 x 
10"'' rps (c) , 8.3 X 10"^ rps (d) and 8.3 x 10"^ rps (e). (f) 
The Newtonian flow field aj]v(r, z) calculated with the finite 
element method (see text). Note the similarity of (a) to (b) 
and (e) to (f). (g) Scatter plot comparison of iooir, z) (a) and 
Los{r, z) for f7 = 8.3 X 10"^ rps (□)(a) and = 8.3 x 10~^ rps 
(o). (h) Scatter plot comparison of ujs{r, z) and the flow field 
of a Newtonian flow uJN{r, z) for fl = 8.3 x 10"^ rps (□) and 
Q, = 8.3 X 10~^ rps (o). (i) vs Q, for comparison to granular 
(A) and Newtonian flow (★). The flow characteristics change 
from granular to Newtonian with increasing shear rate. 



a 



(1) 



Thus, to lowest order, the local stress in a gravitational 
suspension is a linear combination of a frictional stress 
and a purely viscous stress [S] . This is reminiscent of the 
rheology of a Bingham fluid, in that slow flows are rate 
independent while faster flows become dominated by sim- 
ple viscous drag. There is, however, a crucial difference: 
for slow driving, the shear stresses are predicted to be 
proportional to the pressure, while only for faster flows, 
the shear stresses become asymptotically independent of 
pressure. 

Flow profiles — In Fig. [2^-f we compare, for a range 



of driving rates, the measured flow fields, a;s(r, z) (panel 
b-e) with predicted flow fields for dry granular media 
(panel a) and Newtonian flow (panel f). We flx the par- 
ticle flUing height at 23 mm {H/Rs ~ 0.5). Clearly the 
flow structure progressively changes for faster flow rates, 
as the defining characteristic of slow fiows, the trumpet- 
like co-rotating inner core, disappears completely. This 
change is qualitative in nature, with a transition from 
concave to convex shapes of the iso-velocity lines. In ad- 
dition we note an increase of slip near the driving disk 
— while for slow fiows, the normalized angular velocity 
u!s reaches 1 near the disk, for the fastest fiow us has a 
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FIG. 3: (Color online) (a) Rheology of 4.6 mm acrylic par- 
ticles in pure Triton X-100. T{Q,) for a H/Rs = 0.24, 0.51, 
0.60, 0.67, 0.82, 0.91, 1.0, 1.1, 1.2; color (intensity) indicates 
H/Rs. n = 3.0 X 10"* to 0.3. The curve is a fit of the 
form T = To + CSl (see text). The four arrows indicate the 
driving rates where flow profiles were measured, (b) To{H), 
the plateau values as a function H compared to the predic- 
tion from Eq. |2] with /i = 0.57 ± 0.03 for the dry (o) and 
H = 0.59 ± 0.03 for the suspension (-|-) case. 



maximum of 0.7. 

The predicted flow field uJoif, z) for slow dry flows with 
H/ Rs 0.5 is shown in Fig. [2^ - see equations 1,2,6 and 
7 from Ref 25 . The similarity to the slowest flow profile. 
Fig. [2}d, = 8.3 X 10~^ rps, is striking, and is confirmed 
in a scatter plot of a;5(r, z) vs a;£)(r, z), where all data 
for Vl = 8.3 X 10~^ rps (square) collapses on a straight 
line — see Fig. [2^. We conclude that the flow proflles of 
slowly sheared gravitational suspension and dry granular 
media are indistinguishable. 

The predicted flow fleld CLijv(r, z) for Newtonian flows 
with H / Rs ~ 0.5 is determined by a finite element 
software package (COMSOL) to solve the steady state 
Navier-Stokes equations for an incompressible Newtonian 
fluid, and is shown in Fig. [2^. The similarity between the 
measured suspension flows for large $7 and the Newtonian 
flow is, again, striking, and is conflrmed in a scatter plot 
of Ljjs{r, z) vs uJN{r, z), where all data for 51 = 8.3 x 10^^ 
rps (circles) tends to a straight line (Fig. [2ji). 

The crossover from frictional to viscous behavior can 
be quantified further by calculating as function of Q the 
total mean squared deviation (x^) obtained from a linear 
fit of the measured flow profiles ojs{r, z) to the predicted 
dry {u!Dir,z)) and viscous (wAr(r, z)) flows, as shown in 
Fig. [2}. We conclude that the flow profiles of gravita- 
tional suspensions show a crossover from frictional, gran- 
ular behavior to viscous flow upon increasing the driving 
rate. 

Rheology — Can we find the same crossover between 
these two regimes in the rheology? We measure the av- 
erage driving torque as a function of filling height H and 
driving rate il, in order to connect the rheology to the 
findings for the flow profiles discussed above j27 . Since 
index matching is not necessary, we use pure Triton X- 
100 as interstitial fluid; we use the same particles as be- 
fore and keep the temperature fixed at 25°. In each ex- 



periment, ft is incremented from low to high values; each 
data point is obtained by averaging over three or more 
rotations (transients occur over much smaller strains). 

In Fig. [3^ we show T{H, fl) for several different suspen- 
sion filling heights. We conclude that the trends in the 
rheology are similar for all filling heights. First, we ob- 
serve a rate independent regime at small 51, which corre- 
sponds to the range where we observed flow proflles sim- 
ilar to the dry case. Moreover, the overall stress depends 
on filling height, which we will show below to be con- 
sistent with a pressure dependence. Second, the stresses 
become rate dependent for Jl ~ 0.01 rps, and for larger 
rotation rates, the torque increases linearly with 51; over 
the whole range of driving rates explored, the rheological 
data can be well fitted as T = Tq + Cfl, consistent with 
the behavior predicted by Eq. ([T]) _27J. We note here that 
a comparison of the measured torque for pure Triton and 
for the suspension yields that the effective viscosity of the 
suspension is only three to five times larger than r]f. This 
is far below than what would naively be expected from 
textbook formulae, e.g., Krieger-Dougherty. We have no 
explanation for this, but note that in the nontrivial split 
bottom geometry, the suspension packing fraction varies 
throughout the material [25] , which complicates the anal- 
ysis. 

We will now show that the height dependent torque 
for slow flows, To{H), is well described by a prediction 
originally developed for slow dry flows (Fig. [3}d). From 
Eq[l] it follows that the rheology should be determined 
by the local hydrostatic pressure and an effective friction 
coefficient fiQ. Unger and coworkers [33] used these in- 
gredients to predict r(z), the center of the shear band of 
the dry split-bottom flow proflles, but their model also 
gives a prediction for T{H): 

nH 

T{H) = 2gTrp(f>no / {H - z)rVl + {dr/dz)^dz. (2) 
Jq 

Here p is the density of the particles, corrected for buoy- 
ancy in case of submersed particles, (j) is the average pack- 
ing fraction (^ 0.59 [T^]) and /io is the effective friction 
coefficient. Minimization of Eq. [2] yields a prediction for 
T{H) which has not been tested previously. 

As shown in Fig. [Sjs this prediction agrees very well 
with our measurements. The single fit parameter in the 
model allows to accurately extract a friction coefficient, 
which we estimate as po ~ 0.59 ± 0.03. We carried out 
the same measurement of T{H) on dry acrylic particles 
(Fig. [sjs) and obtain a friction coefficient of /Lto = 0.57 ± 
0.03. The two friction coefficients are identical to within 
the experimental error, a fact also observed in Ref. |17j . 
This is strong evidence that in the slow driving rate limit 
the suspension behaves as a dry granular material and 
that lubrication and other hydrodynamic effects can be 
ignored. Furthermore, we can conclude that the simple 
frictional model by Unger correctly captures the overall 
stresses. 
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FIG. 4: (Color online) The rheology of glass beads in a glyc- 
erol mixture at different temperatures. The inset shows the 
same data, with the abscissa rescaled with the viscosity of the 
glycerol at the given temperatures. 



We have also tested the scaling in the viscous regime, 
by measuring the rheology of glass beads [p = 2.5 xlO'^ 
kg/m'^) immersed to H/Rs = 0.4 in glycerol for tempera- 
tures between 4 and 37 °C. The viscosity of the glycerol 
mixture varies more than a decade over this tempera- 
ture range, and hence should change the rotation rate at 
which the viscous regime sets in. The results are shown 
in Fig. |4j Eq. ([T]) requires that the data can be rescaled 
with the viscosity of the liquid rjf — this is indeed ob- 
served in the inset of Fig. [4] Note that the growth of 
torque with strain rate over the larger range probed here 
is somewhat slower than the simple linear prediction |27j . 

Conclusions — Our main finding is that with increas- 
ing shear rate, a gravitational suspension crosses over 
from flowing like a dry granular material to flowing like 
a viscous liquid, consistent with recent modelling of sus- 
pensions based on the inertial number approach. We 
observe this both in the full three dimensional flow pro- 
file, which we revealed using an index matched scanning 
technique, and in rheological measurements. Most of our 
data can be understood based on simple scaling argu- 
ments (to obtain the "transition" shear rate) or elegant 
minimization principles (to obtain p, from T{H)). Our 
measurements indicate that the shape and width of the 
shear band in slow suspensions are the same as for slow 
dry granular flows. Whatever the physics beyond fric- 
tion necessary to produce these flow profiles, our data 
shows that it is equally present in both dry granular and 
gravitational suspension flows. Still, a simple physical 
argument for the most prominent feature of split-bottom 
shear flows — the large width and error function shape 
of the shear zone [T9l - f23] — remains elusive. 
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